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Abstract 

In this paper we describe a search for four radiative decay modes of the D° 
meson: D° -» <frf, D° -> uj, D° -» jf*7, and £>° -> p°j. We obtain 90% CL 
upper limits on the branching ratios of these modes of 1.9 x 10 -4 , 2.4 x 10 -4 , 
7.6 x 10~ 4 and 2.4 x 10~ 4 respectively. 
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I. INTRODUCTION 



Motivated by the successful CLEO II search for b — ► 57 decays Q, we have looked for 
analogous decays in the charm sector. In this paper we consider decays of the pseudo-scalar 
D° meson to final states consisting of a vector meson (<j),u;,K*or p) plus a photon. 

Unlike h — > sj decays, the short-range amplitudes relevant to c — > wy are expected to be 
overwhelmed by much larger long-range electro-magnetic effects. The dominant diagrams 
describing these electro-magnetic amplitudes are shown in Figure [I]. In each case, a pair 
of vector mesons is produced. Providing the quantum numbers are correct, one of these 
can couple to a photon. The phenomenology of such interactions, called "Vector Meson 
Dominance" (VMD), has been well studied ||. Using VMD, one can make rough estimates 
of the expected rates for the modes studied in this paper. If the coupling of the photon to the 
transverse component of a p° results in a vector conversion with about 1% probability, we can 
use the Particle Data Group || value for the D° — > <fip branching ratio, (2.6 ± 0.8) x 10 -3 , 
and expect that BR(D° — > is about 2.6 x 10 -5 • fa, where fx is the fraction of p's 
produced in the decay of the D° which are transversely polarized. Detailed calculations of 
the long-range, lU-exchange and other contributing processes have been published by several 
groups 0-§. The predictions range from 10 -4 to 10~ 6 and are listed in Table |. 

In the b sector, observation of the decay B — > K*^ at the measured rate provided 
compelling evidence for the existence of a "penguin" contribution to the B mesons decay 
amplitude into this channel. The analogous short range penguin diagrams for the radiative 
decay of D° mesons are expected to contribute at the level of S C _> U7 = 1CT 11 — 1CT 8 [J§,(|, 
making them relatively unimportant. 

The long range electro-magnetic contributions that are expected to dominate D° —>■ V7 
decay amplitudes also contribute in the b sector. Their contribution to B — > K*j, for 



example, may be as big as 20% [TC]. It is hoped that a study of these effects in the charm 



sector can improve our understanding of their relevance to bottom decay. 

The CLEO collaboration has recently published a complimentary analysis searching for 
flavor changing neutral currents in D° — > Xl + l~ decays [O . 



II. DATASET AND EVENT SELECTION 

The data used for the analysis described in this paper were acquired with the CLEO II 
detector Jl2| at the Cornell Electron Storage Ring (CESR), and represents a total integrated 
luminosity of 4.8fb~ 1 . 

When searching for D° — > V7 decays we apply several selection criteria on both the 
photon and vector meson candidates before attempting to reconstruct the D° We look 



for — > K + K~ and require 1010 < M^(MeV/c 2 ) < 1030. We also demand that the time of 
flight and specific ionization of both <j) daughter tracks be consistent with Kaon hypotheses. 
We require to candidates to decay into 7r + 7r~7r° and have 763 < M 37r (MeV/c 2 ) < 801. Both 
photons from the 7r° are required to be in the central region of the detector, |cos(0 7 )| < 0.71, 
and the 77 invariant mass must be consistent with a 7r° with x 2 < 4.8. To improve the 
measurement of the 7r° 4-vector, the photons are kinematically fit to the known 7r° mass. 
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We look for K*° -> K-tt + and require 842 < M KlT (MeV/c 2 ) < 942. In this mode we 
also make a cut on the decay angle of the daughter particles in the K* rest frame, requiring 
|cos(#u 7 )| < 0.8, since signal events should follow a sin 2 (6*„ 7 ) distribution due to angular 
momentum conservation. Finally, we reconstruct p's through the decay p — > 7r + 7r~ and 
require 620 < M nw (MeV/c 2 ) < 920. 

In all cases we require that the "radiative" photon be in the central region of the calorime- 
ter, have an energy greater than 830 MeV, and have a calorimeter shower isolated from 
charged tracks in the event. To avoid background from ir° decays we veto photons that are 
part of a n° candidate with x 2 < 15.3. 

In this analysis, all D° candidates are required to come from a D* + — > D°tt + decay. 
The additional kinematic constraint provided by the D* is used to significantly reduce the 
otherwise large combinatoric background. We require the reconstructed mass difference 
between the D* + and the D°, AM = M(D*+) - M(D°), to be between 144.3 MeV/c 2 and 
146.5 MeV/c 2 . To further reduce the background, we demand that Xp* > 0.625, where Xd* 
is defined to be the momentum of the candidate D* divided by the maximum possible D* 
momentum. 

The specific values of the cuts discussed above were chosen after a systematic study of 
{signal) 2 / '(background) for each of the modes, using large samples of GEANT [|l4j] based 
Monte-Carlo data to model each specific signal as well as the background. 



III. BACKGROUNDS 

To learn about possible sources of background for each of the four decay modes, a large 
sample of Monte-Carlo generated e + e~ — > qq events was analyzed. The predominant back- 
ground source found was real D* + — > D°tt + decays where the D° decayed in channels in- 
volving 7r°'s, which in turn decayed such that one of the photons had very little energy and 
went undetected. Since the D* decay in the above sequence is real, backgrounds of this kind 
will result in a false signal that peaks in the mass difference (AM) distribution. Additional 
peaking in the D° mass spectrum will depend on kinematics. 

This type of background is most severe for the D° — > K*^ analysis because poorly 
reconstructed D° — > K~7t + tt° decays, where one of the 7r°'s photons is missed, will peak in 
the D° signal region. Figure shows the D° mass distribution for a set of D° — > K^n + n° 
events analyzed as D° — > K*^. 

In the case of D° — > p7 the problem is less severe since there is no background decay mode 
which peaks in the signal region of our invariant mass distribution, although misreconstructed 
D° — > K~7T + 7r° events cause the upward distortion of the D° invariant mass spectrum just 
below the expected D° mass. Figure |3| shows the distribution of these events when analyzed 
as D° — > p'-f. 

For the modes D° — > wy and D° — > 07 there are no D° decay modes with large enough 
branching ratios to cause noticeable peaking in the reconstructed invariant mass distribution, 
hence we expect the background in the D° mass spectra of these to be smooth. 
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IV. SIGNAL YIELDS AND LIMITS 



All yields were obtained by fitting the D° mass spectra. The signal in all cases was 
parameterized by a double bifurcated Gaussian whose mean and width were determined 
using Monte-Carlo. The background shape used depended on the mode. In the cases of 
D° — > 07 and D° — > wy the background is expected to be smooth and likelihood fits 
were done using simple linear background. The data and fits for these modes are shown in 
Figures | and 

In the cases of D° — > p7 and D° — > K*^, we know the background shape is significantly 
modified by misreconstructed D° — > K^tt + 7t° decays. Using Monte-Carlo, we determined 
the magnitude and shape of this contribution to the D° invariant mass spectrum, and in both 
cases included an additional component in our fits to compensate. The absolute normaliza- 
tion of this additional component was determined from a previous analysis of D° — > K~7t + 7t 
decays |j| . Figures || and |7|(b) show the mass spectra and fits for these modes after sub- 



tracting the contribution from misreconstructed D° —>■ K^tt + ti . 

The results are summarized in Table |. The efficiency for each mode was determined by 
analyzing samples of GEANT based Monte-Carlo "signal" events, and is also presented 
in Table |. To obtain branching ratios from the efficiency corrected yields we performed a 
parallel analysis looking for D* + — > D°tt + , D° — > K~n + decays. Our yield in this mode was 
13, 077 ± 124 events with an overall analysis efficiency of (16.9 ± 0.2) %, determined using 
Monte-Carlo. Using the PDG value of (3.86 ± 0.14) % for the D° — > K~tc + branching ratio 
we find the initial number of D* + — > D°tt + decays in our data sample was (2005 ± 77) x 10 3 . 



V. SYSTEMATIC ERRORS 

Several sources of possible systematic error were investigated, and the results are pre- 
sented in Table ||. With the exception of D° — > wy, the uncertainty in each case is domi- 
nated by uncertainties in fitting. To investigate this error we systematically changed either 
the combinatorial background shape, the normalization of the D° — > K~tt + 7t background 
component (in the pq and K*^ cases only), the signal shape, and the number of bins used 
in the fits. Constant, linear and quadratic background functions were tried. Signal shapes 
were parameterized by Gaussian, Double Gaussian, bifurcated Gaussian and the double bi- 
furcated Gaussian shapes. In each case we took the largest variation as our estimate of the 
systematic error. 

As an additional check we excluded the signal region and fit only the background, using 
simple event counting in the signal region combined with Poisson statistics to obtain the 
upper limits. The result of this procedure for D° — > K* / -f, the mode having the otherwise 
biggest fitting uncertainties, is shown in Figure 0(a). In this case we fitted the D° — > 
K* / -f data with a linear combinatorial component plus the absolutely normalized Monte- 
Carlo predicted D° — > K^7t + tt background, excluding the region between 1.75 GeV/c 2 and 
1.90 GeV/c 2 from the fit. We then count data and predicted background events in the same 
region to obtain a net yield of —33 ± 24. Using a conservative yield of ± 24 events results 
in a 90% CL upper limit yield of 39 events, consistent with the original fitted result. 
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The vector meson mass cuts were studied by varying them to produce a 10% change in 
efficiency and reanalyzing both data and Monte-Carlo with the new values to estimate the 
systematic error. 

To estimate the errors associated with analysis requirements common to all of the studied 
modes (the D* — D mass difference and D* scaled momentum) while avoiding the problem of 
low statistics in the modes of interest, we used numbers obtained in a previous measurement 
of D° -> K-it + it 0]. 

From a CLEO study of the decays rj — > 77 and rj — > 7r 7r°7r°, we assign a 5.5% systematic 
error for uncertainty in the overall 7T° finding efficiency and a 2.5% uncertainty for each 
individual photon. The systematic error due to particle identification was estimated by 
removing that cut entirely and noting the change. 

The systematic errors on the yield and the efficiency were treated separately when cal- 
culating the final upper limits. The efficiency and normalization errors were combined in 
quadrature and the efficiency was reduced by the resulting factor. The fitting systematic 
errors were used to increase the 90% CL upper limit yields. 

VI. CONCLUSION 

Using data representing 4.8 fb _1 of integrated luminosity acquired by the CLEO II de- 
tector at the Cornell Electron Storage Ring, we have conducted a search for radiative decay 
modes of the D° meson. The final results for the 90% confidence level upper limit branching 
ratios for the modes studied are: 

B(D° -> 7 ) < 1.9 x 10~ 4 @ 90%CL 

B{D° -> cj 7 ) < 2.4 x 10~ 4 @ 90%CL 
B(D° -> K*i) < 7.6 x 10" 4 @ 90%CL 
B(D° -> pi) < 2.4 x 10^ 4 @ 90%CL 

We note that all of these values are well above the theoretical expectations as shown in 
Table |. We hope that with more data from CESR, KEK and PEP-II B-factories it will be 
possible to provide improved measurements in the future. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. J.P.A., J.R.P., and I.P.J.S. thank the NYI program of 
the NSF, M.S. thanks the PFF program of the NSF, G.E. thanks the Heisenberg Foundation, 
K.K.G., M.S., H.N.N., T.S., and H.Y. thank the OJI program of DOE, J.R.P, K.H., M.S. 
and V.S. thank the A. P. Sloan Foundation, M.S. thanks Research Corporation, and S.D. 
thanks the Swiss National Science Foundation for support. This work was supported by the 
National Science Foundation, the U.S. Department of Energy, and the Natural Sciences and 
Engineering Research Council of Canada. 
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TABLE I. 



The upper limit yields extracted from the likelihood fit and the resulting 90% confidence 
level upper limits on the branching fractions incorporating systematic uncertainties in yield 
and efficiency determination. 



Mode 


D°- 


+ 4n d 


°^u;7 


£)0 


no Qry 


90% CL Upper Limit Yield 
Detection Efficiency (%) 


8.9 7.7 
5.57 ±0.13% 2.10 ±0.05% 


38.5 
5.51 ±0.13% 


21.6 
5.83 ±0.13% 


Branching Fraction 












90% CL Upper Limit 


1.9 x 


10~ 4 2.4 x 10" 4 


7.6 x 10" 4 


2.4 x 10~ 4 


Theoretical Prediction iMl 0. 


01 - 0.34 x 10~ 4 0.01 - 


0.09 x 10" 4 


0.7-8.0 x 10" 4 


0.01 - 0.63 x 10" 






TABLE II. 








Estimated systematic errors 


for the 


four modes. 








Mode 




D° -> 07 


D° -> w 7 


D° -> K*j 


D° -» p7 


Normalization 




3.87% 


3.87% 


O.O ( 70 


Q Q'7'07' 

o.o i/O 


Monte-Carlo Stat. 




2.25% 


2.4% 


2.3% 


2.3% 


Branching Ratio of the Vector 


Meson 


1.2% 


0.8% 


0.1% 


0% 


Photon and tt° Eff. 




2.5% 


8.5% 


2.5% 


2.5% 


Vector-meson Mass Cut 




3.4% 


2.6% 


2.3% 


3.0% 


Other Cuts 




1.5% 


1.5% 


1.5% 


1.5% 


Particle ID 




9.2% 








Yield/Fitting 




11.8% 


7.3% 


38.8% 


23.6% 


Total Systematic Error 




16% 


12% 


39% 


24% 
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FIG. 1. 

Feynman diagrams for the long distance electro-magnetic contributions. 
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FIG. 2. 

The correctly normalized background contribution from D° — > K~n + n° Monte-Carlo events 
to the D° — > K*^j invariant mass distribution (shaded histogram). The solid line shows the 
expected position and shape for real D° — > K*^ events, also determined using Monte-Carlo. 
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FIG. 3. 

The correctly normalized background contribution from D° — > K~n + n° Monte-Carlo events 
to the D° — > p7 invariant mass distribution (shaded histogram). The solid line shows the 
expected position and shape for real D° — > p7 events, also determined using Monte-Carlo. 
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FIG. 4. 



Data and fit for the D° — > 07 decay mode. 
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FIG. 6. 

Data and fit for the D° — > p7 decay mode. This plot shows the data after subtraction of the 
D° — > K~7t + 7t° background estimation from Monte-Carlo. 
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FIG. 7. 

(a) The observed D° — > K*"f data (points with error bars) and Monte-Carlo predicted 
background (solid histogram). 

(b) Data and fit for the D° — > K*^ decay mode after subtraction of the D° — > K~7t + tt° 
background prediction. 
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